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Abstract 

In this paper, we study the three-body recombination rate of a homogeneous dilute Bose gas 
with a Feshbach resonance at zero temperature. The ground state and excitations of this system 
are obtained. The three-body recombination in the ground state is due to the break-up of an 
atom pair in the quantum depletion and the formation of a molecule by an atom from the broken 
pair and an atom from the condensate. The rate of this process is in good agreement with the 
experiment on 23 Na in a wide range of magnetic fields. 
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I. INTRODUCTION 



Two-body and three-body recombinations are the main reasons for particle loss in cold 
Bose gases. Usually the two-body recombination is dominant in the dilute limit and the 
three-body recombination becomes more important when either the particle density or the 
interaction strength increases. In systems with Feshbach resonances, near a resonance the 
scattering length becomes very large. It was found in the experiment on 23 Na !l| that the 
particle-loss rate increases enormously when close to a Feshbach resonance. The three-body 
recombination rate obtained from the experimental data is much higher than an earlier 
theoretical estimate for a simple dilute Bose gas jjj. 

The large particle-loss rate in 23 Na was explained by one of us in terms of the many- 
body effect due to the special ground-state structure of a dilute condensed Bose gas with a 
Feshbach resonance Q|. The ground state contains both the condensate and the quantum 
depletion. One pair of atoms from the quantum depletion can break up. One atom from the 
broken pair and an atom from the condensate can form a molecule. Therefore in the final 
state, these three atoms become a molecule and an atom with opposite momenta. The rate 
of this process was found to be in good agreement with the experimental data close to the 
resonance Q|. 

There have been considerable progress in the study of the three-body recombination in 
few-body systems 0, S, 0, B, 3|. (See Ref. for a complete review on this subject). 
Particularly, in Ref. ja], the large particle-loss rate in 23 Na was explained in terms of 
the few-body physics. One important prediction from these studies is that the three-body 
recombination rate displays periodic behavior due to Efimov states, which has not been 
observed in experiments so far. 

However both theoretical results from Ref. j^] and Ref. j^] can only fit the experimental 
data of the particle-loss rate |l| very close to the resonance. Theoretical explanation is still 
needed for the experimental data away from the resonance. In this paper, we adopt an 
approach which is valid not only close to the resonance, but also quite far away from the 
resonance, and study a dilute homogeneous system at zero temperature. 
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II. THE TWO-CHANNEL MODEL 



9] 



The starting point of our approach is the two-channel model |9| for systems with a 
Feshbach resonance, 

H = — V^ f • W> + — V0 f ■ V0 + — ip^ip 2 - 7(^V 2 + h.c.) + hdtd), (1) 
2m Am 2 

where m is the atom mass, <fi is the molecular-field operator, and if) is the atomic-field 
operator. The magnetic detuning energy is given by h. The coupling constant of atoms in 
absence of the resonance is given by go = A7ih 2 ao/m, where ao is the asymptotic scattering 
length far away from the resonance. The coupling constant between atoms and molecules 
is given by 7, from which we can define a length scale = 27rfr 4 /(m 2 7 2 ), corresponding to 
the effective range in the two-atom scattering problem |lOj ]. 

In the two-channel model, there are conversions between atoms and molecules. It is 
usually convenient to work in the grand canonical ensemble with the grand thermodynamical 
potential given by 

F = #-M^V + 20V), (2) 
where \x is the chemical potential. 

III. THE GROUND STATE STRUCTURE 

In the dilute region, at zero temperature, both the atoms and the molecules condense. 
The condensed atoms and molecules have phase coherence with each other, and together 
form the condensate. The condensate provides the dominant contribution to the grand 
potential, which is given by 

^0 = I l^ol 4 - 7(0o^o + ex.) + h\U 2 ~ MM 2 + 2|0oD, (3) 

where ipo and 0o are the expectation values of the atomic and molecular fields. The value 
of if) Q and 0o are determined from the saddle-point equations 

dF . . , |2 



dijj 

dF, 



{g Q \M- ^0-2700^0 = 0, (4) 
" -(/i-2/i)^-7^ 2 = 0. (5) 
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In the general solution of these saddle-point equations, the molecular field is always in 
phase or out of phase with the square of the atomic field. For convenience, we choose the 
expectation value of the atomic field to be positive, ipo > 0. 

In the region with the repulsive interaction, the detuning energy is negative, h < 0, 
and the chemical potential is positive, \x > 0. The non-trivial solution of the saddle-point 
equations are given by 

< = A (6) 

where g e g = go — 27 2 / (h — 2/x) is the effective coupling constant between atoms. The 
scattering length is proportional to the effective coupling constant, 

a = mg eS / (Airh 2 ) = a (l - - — — ), (8) 

where the width of the resonance is given by A = 2j 2 /go- In the ground state of a dilute 
Bose gas, most of the particles are condensed atoms [2|, ip 2 m n, where n is the total particle 
density. The chemical potential /i is much smaller than other energy quantities which do 
not have density dependence. 

However as in the traditional theory of a dilute Bose gas, the condensate is not the 
complete picture of the ground state. The ground state also contains the quantum depletion 
which is made of pairs of atoms with opposite momenta. To describe the quantum depletion, 
we need to look at the gaussian fluctuation around the condensate. In the long-wavelength 
and low-energy limit, the gaussian fluctuation in the grand potential is given by 

m = Hfa - /* + 2 ^o - jfzf^V* + - 70o)(W>l k + h.c.)l (9) 

where e& is the kinetic energy of the atom, = h 2 k 2 /(2m). The molecular part of the 
fluctuation is not included in equation Q because molecules have finite energy in the long- 
wavelength limit. However through the virtual excitation of a molecule, the atoms acquire 
a diagonal self-energy approximately given by —2^f 2 ipQ/{h — 2/x), which is included in the 
r.h.s. of Eq.©. 

After applying Bogoliubov transformation, we obtain 

8F 2 = C + Y, E A^ (10) 



where the transformation is given by Ck = Ukipu. + Vkip^_ k , with ul = [1 + (e k + (i)/Ek]/2 and 
v% — u\ — 1. The atoms have a gap less phonon mode with the excitation energy given by 



E k = ^e k (e k + 2fi). (11) 

In the ground state \G), there are no phonon excitations, c^\G) = 0, which is only possible 
if in the ground state atoms with opposite wave-vectors k and — k are paired up. The paired 
atoms form the quantum depletion with the depletion energy C given by 

where V is the total volume of the system. As we will discuss in the later sections, the 
ground state obtained here is not the state with the lowest energy, but a metastable state. 
Nonetheless, in the contest of the traditional theory of a dilute Bose gas, we still name this 
state as the ground state in this paper. 

IV. VACUUM RENORMALIZATION TO THE MOLECULAR EXCITATIONS 

Beside the ground state and phonon excitations, the system also contains molecular 
excitations. From equation (0), the bare molecular excitation energy is simply given by the 
detuning energy h in the long- wavelength limit. However, the molecule energy is strongly 
renormalized by the repetitive process of molecules turning into and forming from atoms. 
In the dilute region, the renormalization to the molecule energy due to the many-body 
interaction is much smaller than vacuum renormalization because the gas parameter is much 



smaller than one, y Snna 3 <C 1. In the following, we study the vacuum renormalization in 
detail. 

In the vacuum, the propagator of the molecular field is given by 

G m (q,n) = - r 1 — — — ■ TT, (13) 

il — y — h + 2fi — S m (q, il) + id 

where S m (q, fl) is the self-energy due to the vacuum renormalization. The self-energy dia- 
grams of the molecular propagator are shown in Fig. ^ As discussed in earlier references 



2. 
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the self-energy of molecules is given by 

S m (q, Q) = 2 7 2 J D(q, Q) [1 + g D(q, Q)}, (14) 
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FIG. 1: The self-energy diagrams of the molecular propagator. The solid lines represent the atom 
propagator and the dashed lines represent the molecule propagator. The vertex represents the 
coupling constant go. 



In Ref. |2j, the second term in the r.h.s of equation (|14j) was dropped. As a result, the 
calculations in Ref. |2j can only be applied when the system is close to the resonance. Here 
we keep both terms in the expression of the molecular self-energy in equation (fT^j) . 

As discussed in Ref. j^J , the divergence in the D-function in equation ()15j) can be removed 
following the standard renormalization procedure by introducing a counter term, which 
yields 



where 




(15) 



and the propagator of the atomic field is given by 




(16) 




(17) 



The molecular self-energy is thus given by 




(18) 



where the parameters v and c are given by 



v 




(19) 



c 



8rf 



a. 



"i 



(20) 
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Although the self-energy given by equation (|18p looks rather complicated, the main fea- 
ture of the the molecular propagator is a pole which is given by 

e a [Ju 2 - 4(c + l)h - u} 2 

o = !i_ 2/ i-iv \ _L L (21) 

q 2 p 4(c+l) 2 V ; 

Close to the pole, the propagator can be approximated as 

C m (q,Q) ^ 



where the renormalization factor Z is given by 



1 + c + 



V 



1 + C + 



- 2// - 9 
z/(c + l) 



4(c+ l)/i-z/ 



(22) 



(23) 



In the limit a ^> ao, the molecule energy is approximately given by the shallow bound-state 



e„e, y , * „ as ,vea in Re, Q. In the „ mit that a is ve ry e.ose to a„, the 

molecule energy is approximately given by the detuning energy, Q m h. 



V. THE THREE-BODY RECOMBINATION 



In the repulsive-interaction region, the molecules have lower energy than the atoms. 
Thus the system can always reduce its energy by forming more molecules. The ground 
state obtained so far is not the state with the lowest energy, but a metastable state. The 
three-body recombination is one of the main processes for the system to move away from 
the metastable state. 

In the metastable condensed state, the three-body recombination takes place in the fol- 
lowing manner. One pair of atoms with opposite momenta in the quantum depletion break 
up. One atom from the broken pair can interact with an atom in the condensate to form a 
molecule, leaving the other atom in the broken pair as an excited atom. Thus the final state 
of this process is the ground state plus a molecule and an excited atom. The probability 
density of the three-body recombination process can be obtained by using Fermi's golden 
rule 

2-7T r d^h 

T = -j- J ^|(k|2 7 ^ 4^ k |G)| 2 5(E fc + ^) 
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2T ' (l4k -Al 2 ^vlZ5(E k + Q k 



h J (2vr) 3 

2tt r d 3 k 2 /i 2 3 



» ■/ (2vr) 3 ' 4e 2 

(24) 



16\/37r7 2 Za^ ■r?; 



m\Qo | 3 

where |k) is the final state which is the ground state plus a molecule with wave vector k 
and an excited atom with wave vector — k. 

In the limit that a 3> ao, the three-body recombination rate T is proportional to a 4 j^j], 
with a coefficient about 140 times larger than that in the single channel case [s]. Away 
from this limit, although there is no simple scaling form for the rate T, we can compute 
its numerical value from equation (|24|). Here the theoretical value of the particle- loss rate 
in the 23 Na system is plotted and compared with the experimental data |l| in Fig. El The 
experimental data is the loss rate divided by the product of total particle number N and the 
average of the density squared. The theoretical value is the loss rate divided by the product 
of N and the homogeneous density squared. All the parameters used in the calculation are 
taken from the experiment, a$ = 3.3nm and A = 2\ibG. The theoretical result is fairly 
close to the experimental data in a wide range of magnetic fields in the repulsive-interaction 
region, although the theoretical derivation so far is for the homogeneous system and the 
experimental system is inhomogeneous. 



VI. CONCLUSIONS 

We have studied a homogeneous dilute Bose gas with a Feshbach resonance is studied 
at zero temperature. In the ground state, there are both the condensate and the quantum 
depletion. There are two types of excitations, phonon excitations and molecular excitations. 
The three-body recombination in the ground state is due to the formation of a molecule 
from an atom in the condensate and an atom from the quantum depletion. The excellent 
agreement between our calculation and the experimental results in a wide range of magnetic 
fields leads us to conclude that the many-body effect due to the ground-state structure plays 
a major role in the three-body recombination process. 

An interesting question is whether or not such theoretical description is still valid when 
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FIG. 2: The coefficients of the particle-loss rate for the Na system, N/N(n 2 } is plotted vs the 
magnetic field as in Fig. (2) of Ref. The straight line is the theoretical result computed from 
Eq.pil): the squares and triang 



es are the experimental data taken at ramp speeds 0.13G/ms and 
0.31G/ms of the magnetic field jj. All the parameters used in the calculation are taken from the 



experiment 



D. 



the system moves further away from the resonance. When the system is very far from 
the resonance, the two-channel description should probably be replaced by a single- channel 
model and the three-body recombination rate is likely to recover the single-channel form. 
One possible criterion is the comparison between the molecule energy in the two-channel 
model and the bound-state energy in the single-channel model. If the molecule energy in 
the closed channel is much deeper than the bound-state energy in the open channel, the 
single-channel description is probably better. 
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